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Abstract 

While it is often stated that the notion of electroweak (EW) naturalness in super¬ 
symmetric models is subjective, fuzzy and model-dependent, here we argue the contrary: 
electroweak naturalness can be elevated to a principle which is both objective and pre¬ 
dictive. We demonstrate visually when too much fine-tuning sets in at the electroweak 
scale which corresponds numerically to the measure Nbg ^ ^ 30. While many 

constrained SUSY models are already excluded by this value, we derive updated upper 
bounds on sparticle masses within the two-extra parameter non-universal Higgs model 
(NUHM2). We confirm the classic Barbieri-Giudice (BG) result that Abg < 30 implies 
fi < 350 GeV. However, by combining dependent soft terms which appear as multiples of 
777 . 3/2 in supergravity models, then we obtain mg ^ 4 TeV as opposed to the BG result 

that mg ^ 350 GeV. We compare the NUHM2 results to a similar scan in the pMSSM 
with 19 weak scale parameters. In the pMSSM with complete one-loop scalar potential 
plus dominant two-loop terms, then a mg < 7 TeV bound is found. Our tabulation of 
upper bounds provides a target for experimenters seeking to discover or else falsify the 
existence of weak scale supersymmetry. In an Appendix, we show contributions to the 
naturalness measure from one-loop contributions to the weak scale scalar potential. 
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1 Introduction 


Weak scale supersymmetry [T] (SUSY) is a highly motivated paradigm for physics beyond the 
Standard Model (SM). The principal motivation is that it offers a technical solution to the 
gauge hierarchy problem or naturalness problem of the SM in that it ensures cancellation of 
quadratic divergences of scalar field masses to all orders in perturbation theory [2]. This is 
especially relevant now that a bonafide fundamental scalar has been discovered[3l 0]: i.e. the 
Higgs boson h. The SUSY technical solution is accomplished in a highly simple manner: merely 
extending the set of spacetime symmetries which underlies quantum field theory to their most 
general structure based upon a graded Lie algebra. In fact the SUSY paradigm is supported 
indirectly via three sets of measurements: 


1. the measured values of the three gauge forces at the weak scale are exactly what is needed 
for SUSY gauge coupling unification [5], 

2. the measured value of the top quark mass is exactly what is needed to properly drive a 
radiative breakdown in electroweak symmetry [H] and 

3. the measured value of the Higgs boson mass mh = 125.09 ±0.24 GeV falls squarely within 
the predicted narrow band of allowed MSSM values [7] . 


On the negative side of the ledger, there is as of yet no sign of supersymmetric matter 
after extensive runs of LHC at ^/s = 7 — 8 TeV[8l [9]. This fact has led many theorists to 
question the validity of SUSY in light of the oft-repeated mantra that weak scale naturalness 
requires weak scale sparticles. It is also baffling on the experimental side as to when the weak 
scale SUSY hypothesis is ruled out, and when one ought to move on to alternative directions. 
This is especially vexing in that many theoretical predictions tend to he just beyond current 
exclusion limits: when the exclusion limits increase, then the theoretical predictions retreat 
towards higher mass values again just beyond the latest lower mass bounds. This leads to the 
important question: what does it take to falsify the weak scale SUSY hypothesis! When is the 
job done? Are new accelerators and experiments required, or does LHC with ^/s ~ 13 — 14 
TeV and high luminosity have the necessary resolving power? 

Historically, an answer to this question was provided in the classic paper by Barbieri and 
Giudice[T0] wherein- based upon the “naturalness criterion”- they derived upper bounds on 
various sparticle masses Q These upper bounds could serve as targets for experimental facilities 
with the intent to discover or disprove weak scale SUSY. Following earlier work by Ellis et 
a/. [12], they introduced the naturalness measur^ 


Abg = maXi 


d log m| 
d log Pi 


( 1 ) 


where the Pi are fundamental parameters of the theory labelled by the index i. Working within 
the MSSM with unified GUT scale soft breaking terms mo, mi/ 2 ,Ao and B, they presented 

^See also Ref. m- 

^The authors of Ref’s [13] note the link between naturalness and Bayesian statistics wherein naturalness 
corresponds to “more probable”. 
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upper bounds on sparticle masses and parameters as a function of the top-quark mass (which 
of course was not yet known at the time). In light of our present knowledge of the top quark 
mass {rrit = 173.21 ± 0.87 GeV PDG value with combined statistical/systematic errors), they 
would conclude that in order to accommodate Abg < 10 (or better than 10% fine-tuning in 
mz), then mo ^ 300 GeV, mi /2 ~ 100 GeV and /i ~ 200 GeV. These upper bounds implied 
that the lighter charginos m^^ ~ 100 GeV and gluinos rrig ~ 350 GeV. Explorations at LEP2 
resulted in limits of m^^ > 103.5 GeV[T6] so that weak scale SUSY already appeared somewhat 
fine-tuned by the post-LEP and pre-LHC era[IZ|. The current limits on gluino mass from LHC8 
require rrig ~ 1300 GeV. Thus, in the post LHC8 era, the naturalness issue has intensified, and 
based upon these theory/experiment confrontations, one might well be tempted to conclude 
that weak scale SUSY has been disproved [IHj. 

In this paper, we re-examine upper bounds on sparticle masses from the naturalness princi¬ 
ple. In Sec. we discuss naturalness and fine-tuning and articulate the Naturalness Principle. 
In Sec. 1^ we apply the Naturalness Principle at the weak scale to derive the electroweak fine- 
tuning measure Aew- In Sec. |^we discuss Barbieri-Giudice fine-tuning Abg, how it depends 
on the selection of an independent parameter set and how it relates to Aew- when properly 
applied by combining dependent contributions, then Abg — Aew- In Sec. we demonstrate 
visually when too much fine-tuning sets in and what values of Aew — Abg are too much. In 
Sec. [^we derive upper bounds on sparticle masses in the NUHM2 model and compare them to 
upper bounds from the classic BG paper. In Sec. we derive alternative upper bounds arising 
from a scan over the 19 dimensional weak scale pMSSM parameter space. We compare these 
against results from the 19 parameter SUGRA model with soft terms defined at Q = tugut- 
In Sec. 1^ we conclude by presenting a bar chart of upper bounds on sparticle masses which 
form a target for experimenters seeking to confirm or refute weak scale supersymmetry. The 
hard target is that /i < 350 GeV in any models (defined at a high scale or weak scale) based 
on the MSSM. In models with RG running from mauT to rriyjeak, then rrig ~ 2 (4) TeV for 
Abg — Aew < 10 (30). If one dispenses with RG running as in the pMSSM, then no bound 
is obtained on rrig using the one-loop RG-improved effective potential. By including leading 
two-loop contrubutions to the scalar potential, we find a bound of rrig < 7 TeV for Aew < 30 
in the pMSSM. In an Appendix, we discuss the contributions to Aew from various radiative 
corrections S“(f) and plot their magnitudes in 2-dimensional parameter planes. 

2 The naturalness principle 

Several definitions of naturalness can be found in the literature: some are more abstract while 
others are more pragmatic. Here we articulate the following Naturalness Principl^TI 


^ Dimopoulos and Susskind articulate: “Naturalness: no parameter needs to be adjusted to unreasonable 
accuracy.” m 

^ Weinberg states: “The appearance of fine-tuning in a scientific theory is like a cry of distress from nature, 
complaining that something needs to be better explained” BO]. 


2 



An observable O is natural if all independent contributions to O are less than or 
of order O. 

Suppose O can be calculated in terms of n independent contributions O = Oi + • • • + On- If 
one of the contributions, say o„, is far larger than O, then it would have to be the case that 
one or more of the o* would have to be a large opposite sign contribution which would require 
fine-tuning of o* ~ —o„ such that 0 <C o„. This is the link between naturalness and fine- 
tuning: a quantity O is natural if it requires no large fine-tuning of independent contributions 
to maintain its measured value. 

A common pitfall in evaluating when a quantity is natural is to split it into dependent 
contributions: if 0 = 0i -|- 02 -|- • • •, but as one increases oi then 02 necessarily increases to large 
opposite-sign, then oi and 02 should be combined before the final evaluation of naturalness. 
This pitfall has been dubbed the “fine-tuning rule” BB 

To see how is a measure of naturalness, consider the observable O expressed as a linear 
combination of n fundamental parameters pp 


O — aipi -!-••• + anPn ( 2 ) 

where the a* are numerical co-efficients. Applying the Abg measure, we would find 

= maXi \aiPi/0\ (3) 

i.e. the BG measure just picks off each term on the right-hand-side of Eq. and compares 
it to To avoid fine-tunings, then Abg should be less than some value (typically 10-100) 
depending on how much fine-tuning one is willing to tolerate. 


Abg = rnaxi 


Pi do 


Odp, 


3 Weak scale naturalness 


Starting with the weak scale scalar (Higgs) potential of the MSSM, the minimization conditions 
= 0 allow one to determine the Higgs field VEVs in terms of the soft SUSY breaking 

parameters and the p parameter[T]. Then, since m| = + v1)/2, we can relate the 

observed value of to the weak scale SUSY parameters as 


ml 

2 


+ ^d) - + ^u) tan^ ^ 

(tan^ /3 — 1) 


~ -m% - 


(4) 


Here, and are the weak scale soft SUSY breaking Higgs masses, /i is the supersym¬ 
metric higgsino mass term and S“ and contain an assortment of loop corrections to the 
effective potential]^ Already at this stage (not yet worrying about high scale parameters), the 
naturalness principle requires each term on the right-hand-side (RHS) of Eq. |^to be compara¬ 
ble to or less than m|/2. Erom the partial equality on the RHS, it is plain to see that the weak 

®If terms include powers of Pi, e.g. O = then each sensitivity co-efficient contains an additional 

factor of Uj. 

®We will not consider tan /3 as an independent parameter here. Its value is determined also by the minimiza¬ 
tion conditions in terms of the SUSY parameters. 
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scale values of rri\j , fj? and the various S“(i) {i labels the various loop contributions) should 
all be comparable to or less than m|/2[22]. This allows us to define the electroweak fine-tuning 
measure /S.ew as 


^Ew = max \each term on RHS of Eq. /(m|/2). (5) 

In gravity-mediation, the GUT scale soft terms are expected to be of order the gravitino 
mass m 3 / 2 . Since the physical sparticle masses are derived from the soft terms, we thus expect 
the weak scale soft terms also to be of order 771 . 3 / 2 . Since LHC8 requires multi-TeV values of 
mg and m^, then it seems LHC8 is telling us that 7773/2 is also in the multi-TeV range. This is 
good news for the SUSY flavor and CP problems [23] and the gravitino problem |24j: all these 
are ameliorated by multi-TeV values of 7773/2 and sparticle masses. The puzzle then is: why is 
mz (and mw and mif) not also at the multi-TeV scale? There is one soft term whose weak 
scale value may be very different from its GUT scale value. The term m^ may hold multi-TeV 
values at mcuT but is necessarily driven through zero to negative values at Q = m^eak in order 
to break electroweak symmetry [6] . If it is driven to small values comparable to —m‘zl2 rather 
than large negative values, then one naturalness condition may be satisfied. This case has been 
dubbed radiatively-driven na,tura,lness fIf\YIE\ . 

The other condition for naturalness is that iJ? ~ 777|/2|^ In this case, we note that /7 is SUSY 
conserving and not SUSY breaking. It is a very different entity from the soft terms: naively, it 
would be present even in the absence of SUSY breaking. In this case, one expects its value to be 
of order the GUT or Planck scale Mp. How it comes instead to be ~ mz is known as the SUSY 
/7 problem [TT] . There are two parts to the solution of the SUSY /r problem: first, one must 
forbid it (via some symmetry) from attaining values ~ Mp, and second, one must regenerate 
it of order mz (via symmetry breaking). In the original Kim-Nilles formulation[SI|, it was 
noted that in the SUSY DFSZ[32| axion model, Peccei-Quinn symmetry forbids the /7 term. 
The spontaneous breaking of PQ symmetry at scale /„ ~ 10^^ GeV generates an axion (thus 
solving the notorious strong CP problem), but also generates a mu term of order /a/Mp^Sj- 
Since one expects the gravitino mass (and hence soft terms) to be of order 7773/2 ~ ^‘hidden/ 
where mhidden is an intermediate mass scale associated with hidden sector SUGRA breaking, 
then the apparent Little Hierarchy fi <C 7773/2 is just a consequence of a mismatch between 
PQ breaking scale and hidden sector mass scale: fa mhidden- In fact, in models such as the 
Murayama-Suzuki-Yanagida (MSY) SUSY axion model[33) PQ symmetry is radiatively broken 
as a consequence of SUSY breaking. In the MSY model, for canonical parameter values a small 
value of /i ~ 100 — 200 GeV can be easily generated from 7773/2 "^5 — 20 TeVjHSj. 

Finally, the radiative corrections S“ should be less than or ~ m^/2. Typically, the top 
squark contributions are the largest of these. The top squark contributions S“(ti 2 ) are min¬ 
imized for TeV-scale highly mixed top squarks, which also lift the Higgs mass to mh ~ 125 

GeV [2S|. 

The EW fine-tuning measure is the most conservative of the fine-tuning measures in that 
any model with large Apw is surely fine-tuned. Also, Apw clearly agrees with Apo for SUSY 

^In this regard, we rely on Einstein’s advice to maintain the theory as simple as possible, but no simpler. 
By adding various exotica to the MSSM, then one can create models with heavy higgsinos which may still be 
natural [271 UHl US] . These often involve adding exotic states such as scalar gluons. It is not clear whether such 
constructs admit a UV-completion. 
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models defined purely at the weak scale such as the pMSSM. Another virtue of IS.ew 

is that its value is model-independent (within the MSSM) in that it doesn’t depend on how 
the weak scale spectrum was generated. Also, from a pragmatic point of view, the fine-tuning 
encapsulated in Aew is exactly where spectrum generators invoke explicit fine-tuning. In such 
codes, usually and other soft terms are calculated at the weak scale via RG running, and 
then the value of /i is dialed/fine-tuned to enforce that the measured value of is obtained. 
Without such fine-tuning, the generated value of mz would typically lie in the multi-TeV 
region| 2 T]. 

But does Aew encapsulate all the fine-tuning, including high-scale effects, or is it just a 
lower bound on fine-tuning[3H]? 

4 BG fine-tuning and independent model parameters 

For models defined in terms of high scale parameters, the BG measure can be evaluated by 
expanding the terms on the RHS of Eq. using semi-analytic RG solutions in terms of fun¬ 
damental high scale parameters uni EH]. For the case of tan /? = 10 and taking the high scale 
A = mQUT, then one finds [101SH 02] 

m| ~ -2.18^^ + 3.84M| - O.fiSMgAt - 1.27m^^ - 0.053m^^ + 0.73771^3 + 0.57m^3 + • • • (6) 

The problem with most applications of the BG measure is that in any sensible model of 
SUSY breaking, the high scale SUSY parameters are not independent. For instance in gravity- 
mediation, for any given hidden sector, the soft SUSY breaking terms are all calculated as 
numerical co-efficients times the gravitino mass [131 011 03]: e.g. M 3 (A) = aM 3 ?^ 3 / 2 , Aj = 
^At?^ 3 / 2 , 'iTT'Qg = ^Q 3 ^ 3 / 2 ) where the a* are just numerical constants. (For example, in 
string theory with dilaton-dominated SUSY breaking [HUS], then we expect ml — ml 12 with 
mi /2 = — Ao = \/ 3 m 3 / 2 ). The reason one scans multiple SUSY model soft term parameters is 
to account for a wide variety of hidden sectors possibilities. But this doesn’t mean each soft 
term is independent from the others. By writing the soft terms in Eq. [^as suitable multiples of 
ml 12 , then large positive and negative contributions can be combined/cancelled and one arrives 
at the simpler expression [21103] : 


m| ~ -2/r^(A) -F a ■ ml /2 (7) 

since the SUSY mu term hardly evolves. The value of a is just some number which is the sum 
of all the coefficients of the terms oc ml/ 2 . 

Using the BG measure applied to Eq. then it is found that naturalness requires ~ 
/i^(A) ~ m| and also that aml /2 ^ The first requirement is the same as in Aew- The 
second requirement is fulfilled either by mj ,/2 ^ mz (which seems unlikely in light of LHG 
Higgs mass measurement and sparticle mass bounds) or when 7713/2 is large but the co-efficient 
o is small 1211 : te. there are large cancellations in Eq. By equating in terms of weak 

^The pMSSM, or phenomenological MSSM, is the MSSM defined with weak scale input parameters where 
all CP violating and flavor violating soft terms have been set to zero. Also, usually first/second generation soft 
terms are set equal to each other to avoid flavor-violations. 
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scale parameters (Eq. Q with m| in terms of high scale parameters Eq. [^, and using the 
fact that /u(A) ~ ii{weak), then also amlj 2 — m\j^{weak) and so a low value of A^g also 
requires a low value of m\^{weak) ~ —m\. By properly evaluating BG fine-tuning in terms 
of independent SUGRA parameters, namely 771 . 3/2 and /i(A), then we are lead back to the same 
sort of conditions as implied by low Aew- he. that is driven radiatively to small negative 
values. In this manner 

• Aew encompasses both high scale and weak scale fine-tuning. 

The ambiguity between fine-tuning measures is removed]^ Eor electroweak naturalness in SUSY 
theories, it is not the case that sparticles need to be near the scale ~ 100 GeV: it is just the 
weak scale Lagrangian parameters and jj,: the remaining soft breaking terms may lie 

comfortably in the multi-TeV range at little cost to naturalness. 

5 How much fine-tuning is too much? 

Once a reliable measure of fine-tuning is established, then the next question is: for which values 
of Aew is a model natural, and for which can it be considered fine-tuned? The original BG 
paper considered A^g < 10 to be natural. However, as experimental limits on sparticle masses 
grew, then much higher values of Abg were tolerated- up to ~ 100[l2] or even ~ 1000[50]. This 
increased tolerance perhaps reflected a reluctance to easily give up on an amazingly beautiful 
and simple paradigm even in the face of apparent large fine-tuning. 

By properly evaluating Abg in terms of Eq. or equivalently Aew (which includes one- 
loop radiative corrections), then we can re-evaluate how much fine-tuning is too much. In 
Fig. B we show the dominant contributions to Aew for a simple NUHM2 benchmark model 
with mo = 5 TeV, mi /2 = 700 GeV, Aq — —8.3 TeV and tan/5 = 10 with m^ = 1 TeV and 
variable fi. In the first case, with = 110 GeV and Aew = 13.6, we display the several largest 
contributions to Aew- Some are positive and some are negative but all are comparable to ±10. 
This case is visually highly natural: the Z mass is ~ 100 GeV because its various contributions 
are ~ 100 GeV. In the second column, with /r = 200 GeV and Aew = 22.8, the contribution 
to Aew is largest but is on the whole balanced by several comparable positive terms: again, 
one would not claim it as unnatural. The third column with fi = 300 GeV, we see that the Hu 
contribution to Aew is again the largest, but in this case now a value of fi must be selected as 
large positive to compensate and enforce that mz = 91.2 GeV: this case is starting to already 
become unnatural. By the fourth column with fi = 400 GeV and Aew = 51.5, the fine-tuning 
is visually striking: the model is no longer natural. The unnaturalness is only accentuated in 
columns five and six for /i = 500 and 600 GeV respectively. 

® Another common fine-tuning measure [131311133] is known as Higgs mass or large log fine-tuning Ajjs- In its 
usual implementation, Ahs requires the radiative correction to the Higgs mass ~ 

be comparable to m\. This contribution is usually written as SrrAj^\rad ~ ~^iAAAQ^+'niif^+A^) In {A?/rtigugy) 
which is used to claim that third generation squarks ^ be approximately less than 500 GeV and At be small 
for naturalness. This expression for makes several approximations- the worst of which is to neglect that the 

value of itself contributes to By combining dependent contributions, then instead one requires that 

the two contributions (A) -|- be comparable to m|. Since {A)+8m\^ = rn?H^ [weak), 

then we are lead back to the same conditions as low Aew and low Abg- 
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Figure 1: Plot of contributions to Aew for various values of superpotential parameter. 

To be conservative, it is evident that fine-tuning has set in for fi values ~ 300 — 400 GeV. 
For such cases one would expect the value of mz also to be in the 300 — 400 GeV range. 
With a value of /r = 350 GeV as a conservative upper estimate, and with the contribution 
Aew{iA — /r^/(rrr|/2), we would then expect already values of Aew > 30 to be overly fine- 
tuned. This is somewhat above the values expected in the original BG paper where they 
adopted = 10. If we increase the mass bounds from BG by a factor A^^/IQ then we 

expect from BG with = 173.2 GeV that ji < 350 GeV and rrig ~ 350 GeV. 

From scans over the popular mSUGRA/GMSSM[511 152] model, it is found that with nih = 
125 ± 2 GeV, the lowest value of Aew which can be obtained is ~ 100|53l [SH [2T]. In this 
case, the mSUGRA/GMSSM model is highly fine-tuned and is already ruled out. In addition, 
an assortment of other models- including mGMSB, mAMSB and mirage unification models- 
are also ruled out [21]. What remains is a region of NUHM2 parameter space where small 
fi ~ 100—300 GeV is allowed and where highly-mixed top squarks may live in the few TeV range. 
This region of NUHM2 parameter space is labelled as RNS, standing for models with radiatively- 
driven natural supersymmetry [251126] . RNS models are characterized by light higgsinos with 
mass ~ /i ~ 100 — 300 GeV. The LSP is a higgsino-like WIMP with a thermally-produced 
underabundance of dark matter. However, solving the QCD sector naturalness problem via the 
axion[nH| leads also to axion dark matter[nS| so that one expects two dark matter particles: 
the axion along with a higgsino-like WIMP [57] . Variants of this model with greater parameter 
freedom, such as non-universal gaugino masses leading to natural SUSY with bino-like or wino- 
like LSPs- are also allowed [HR]. 
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6 Mass bounds from naturalness in NUHM2 


To derive sparticle mass bounds from the naturalness principle, we will generate SUSY spectra 
using Isa jet [Sni [03 in the 2-parameter non-universal Higgs model [HI] (NUHM2) which allows 
for very low values of /S.ew < 10. The parameter space is given by 

"^- 0 , ■^ 1 / 2 ) ^ 0 ) tan/3, /i, niA, (NUHM2). (8) 

The NUHM2 spectra and parameter spread versus /S.ew were evaluated in Ref. [SH] but with 
rriA restricted to < 1.5 TeV and mi /2 < 2 TeV. Here, we update these results by ensuring that 
we use sufficiently large range of input parameters that our upper bounds on sparticle masses 
surely come from /S.ew < 30 rather than from artihcial upper limits on scan parameters. Some 
previous mass bounds were extracted in Ref ’s (20111ES]. Here, we enlarge the scan region to 
include: 


mo 

: 0-20 TeV, 

mi/2 

: 0.3-3 TeV, 

-3 

< Ho/mo < 3, 

/i 

: 0.1 - 1.5 TeV, 

rriA 

: 0.15-20 TeV, 

tan P 

: 3-60. 


We require of our solutions that: 

• electroweak symmetry be radiatively broken (REWSB), 

• the nentralino Zx is the lightest MSSM particle, 

• the light chargino mass obeys the model independent LEP2 limit, > 103.5 GeV [T6] . 

• LHC8 search bounds on rrig and niq from the mo vs. mi /2 plane[Hj are respected, 

• m/, = 125 ± 2 GeV. 

It is important to note in this Section that some of our upper bounds come from the specihc 
model we sample. Eor instance, the extracted npper bonnd on the glnino mass comes from RG 
running effects where the gluino mass feeds into the top squark soft terms, and the top-squark 
soft terms are mainly constrained by the E“(ti_ 2 ) contributions to /S.ew- In this respect also, 
/S.EW is sensitive to high scale effects since it measures if a particular model defined by high 
scale parameters can generate the weak scale characterized by mw,z,h ~ 100 GeV. In the next 
Section, we will discuss model-independent upper bounds which do not depend on high scale 
physics. 

The first results of our scan are shown in Eig. In frame a), we plot Aew vs. rrig. The 
symbols are color-coded according to low (< 15), intermediate (15 — 30) and high (> 30) values 
of tan 13. Erom the plot, we see first that there is indeed an upper bound to rrig provided by 
naturalness. If we enforce Aew < 30, then we find that 



• rrig ~ 4 TeV. 

This is to be compared to the reach of LHC for gluinos via gluino pair production followed 
by cascade decays. For the case of heavy squarks rriq rrig and with ~ 1000 fb“^, LHC13 
has a reach up to rrig ~ 2 TeVjHH ES]- Thus, while natural SUSY may well be discovered by 
LHC15p^ there is also plenty of natural SUSY parameter space well beyond the LHC reach for 
gg production. This bound may also be compared to the original BG result for rrit = 173.2 
GeV and with Abg < 30: there it was found that rrig ~ 350 GeV. The discrepancy between 
results arises because BG evaluated Abg within the four-parameter mSUGRA/GMSSM effec¬ 
tive theory where tuq, mi/2, Aq and B are assumed as independent. By recognizing the soft 
terms as dependent multiples of 771 . 3 / 2 , then using Eq. [^we have Abg — Aew and much larger 
values of rrig are allowed while preserving naturalness. Our result may also be compared with 
Feng|l2] who uses an upper bound of 1% hne-tuning in a multi-parameter SUSY effective the¬ 
ory: he then Ends rrig ~ 1.4 TeV, slightly beyond the latest bound of rrig ~ 1.3 TeV from LHC8 
searches. If we scale this back to 3.3% hne-tuning to compare with our result, then Feng would 
obtain rrig ~ 770 GeV, well below the LHG8 lower limit on rrig. 

In Fig. we also show Aew versus b) the higgsino mass /r and c,d) versus the gaugino 
masses Mi and M 2 . Since ~ ^ ~ |/i| in models with a higgsino-like LSP, then we 

expect that 

• 350 GeV 

(the lighter the better). While the higgsino-like electroweak-ino masses are necessarily not too 
far from mz,h, they are notoriously difficult to see at LHG due to their compressed spectrum jHHl 
EH]- A possible way forward may be via pp —?• ZiZ 2 g production followed by Z 2 —?• pA ij,~Z^. In 
this case, the hard gluon ISR serves as a trigger so that events containing soft dimuons may be 
visible above SM backgrounds [66l [671 EH]- In contrast, the required light higgsinos should be 
easily visible at an e+e“ collider operating with ^/s > 2/i[65l [69] where also their masses and 
mixings can be extracted to high precision. In this, we are in accord with the conclusion of 
BGdDI that lepton colliders provide a more powerful probe of SUSY electroweak naturalness 
than do hadron colliders. 

The bound extracted from frame c) is that the bino mass Mi ~ 900 GeV for Aew < 30. 
Since the NUHM 2 model assumes gaugino mass unihcation, this translates to the mass on the 
third lightest neutralino Z 3 which is then mainly bino-like. 

The bound on the wino mass M 2 shown in frame d) translates to a bound on the wino-like 
electroweakinos of rri^^ ~ 1600 GeV. As noted in Ref’s [701 E3| , wino pair production at 

LHG13 (via W^Z/^ and 14^2"14^2“ production) forms the dominant visible reaction. In fact, the 
reach of LHG13 via the same-sign diboson (SSdB) signature (pp —?• 142-^4 4F^4U^+ fJx) 

exceeds the reach via gluino pair production for integrated luminosities r\j 300 fb The reach 
of LHG13 for wino pairs via the SSdB signature extends to ~ 680 GeV. 

In Fig. 1^ we show updated upper bounds on third generation squark masses. From frame 
a), we see that 

• mp ~ 3 TeV. 

^*^For an overview of what natural SUSY looks like at LHC13, see Ref. |S3| 
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Figure 2: Plot of /S.ew vs. rrig, /i, M\ and M 2 from a scan over the NUHM2 model. Points 
with /S.EW < 30 (below dotted line) are considered natural. 
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Figure 3: Plot of ^ew vs. mi^ and from a scan over the NUHM2 model. 


This upper bound is much higher than previous incarnations of natural SUSY based on large-log 
fine-tuning [inj (where it was claimed that naturalness required three third generation squarks 
with mass ~ 600 GeV). The projected reach of LHC13 for rrif^ in various simplified 

models extends up to mq ~ 1 TeV. Thus, the lighter top squark may also lie well beyond 
LHC13 search capabilities with little cost to natnralness. From frames b), c) and d), we find 
successively that rrif^ ~ 9 TeV, ~ 9 TeV and ~ 10 TeV. 

In Fig. lilwe show plots of Aew versus various matter scalar masses a) b) c) nif^ 
and d) m^ilFrame a) is typical of all first/second generation matter scalars. Here, we extract 
an upper bound of ~ 10 TeV. The bound on first/second generation matter scalars arises 
from D-term contributions to S“[26]. For certain mass degeneracy patterns listed in Ref. ini, 
these contribntions nearly cancel amongst themselves. The limits become much stronger for 
non-degenerate matter scalars. The limits are also affected by two-loop RG effects where heavy 
first/second generation matter scalars feed into third generation and Higgs soft term evolution. 
In this case, large mo(l, 2) can drive to large instead of small negative values [7H]. The limit 
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Figure 4: Plot of Aew vs. rriu^, rufj and from a scan over the NUHM2 model. 


on first/second generation sleptons, shown in frame b), is similar. If we allow for non-degenerate 
generations, i.e. mo(l,2) 7 ^ mo(3), then the upper bounds on first/second generation squarks 
and sleptons can increase to ~ 20 TeV|2S]- 

In frame c), we show Aew vs. mfj. In this case, the upper bound on third generation 
sleptons is also nif^ ~ 10 TeV. An npper bound on ttia can be extracted from frame d) of Fig. 
[4 Here we find ttia ~ 5 (8) TeV for tan/3 < 15 (50). Thus, the heavy Higgs bosons with mass 
niH ~ m^f± ~ niA may also be well beyond the reach of LHC13 at little cost to naturalness. 

6.1 Comparison of original BG results to this work 

In Table [T} we compare the upper limits on sparticle masses and the /r parameter extracted 
from the original BG paper[T0] as compared to this work. The BG results are presented for 
rrit = 173.2 GeV and scaled to impose Abg < 30 instead of Abg < 10. The BG results used 
mo, mi/ 2 , Aq and fi as independent parameters whereas our results combine these contributions 
to m| since each is computed as a multiple of ms /2 in gravity mediation models. To include 
radiative corrections, we use the Aew measure. 
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Since both groups use /i as an independent parameter, then both groups agree on the upper 
bound jjL ~ 350 GeV leading to relatively light higgsinos: the closer to mh the better. However, 
for other sparticle masses, then our results differ markedly. Whereas BG find rrig ~ 350 GeV, 
we find nig ~ 4 TeV- possibly well out of range of LHC13. Also, for the bino mass Mi, BG find 
Ml ~ 90 GeV while we find Mi ~ 900 GeV. For the wino mass M 2 , BG find M 2 ~ 170 GeV 
(leading to the suggestion that charginos should appear at LEP 2 if SUSY is not fine-tuned) 
while we find M 2 ~ 1700 GeV (possibly out of range of LHG13 even via wino pair production 
with decays to same-sign dibosons [ZO])- Meanwhile, for matter scalars, BG found ~ 700 
GeV whereas we find rriu^ ~ 10 TeV. If we allow for non-degenerate matter scalar generations 
(mo(l, 2) 7 ^ mo(3)) then our limit can increase to r\j 20 TeVpi)]. Likewise, BG found mg^ ~ 520 
GeV whereas we expect rrie^ ~ 10 TeV (20 TeV for split families). 


mass bound (GeV) 

BG 

this work 

p 

350 

350 

m~g 

350 

4000 

Ml 

90 

900 

M 2 

170 

1700 


700 

10000 ( 20000 ) 

men 

520 

10000 ( 20000 ) 


Table 1 : Upper bounds on masses (in GeV) from naturalness with A < 30 from original BG 
paper compared to Sec. |^of this work in the NUHM2 model. The entries in parenthesis would 
result if one allows for non-degenerate generations of soft scalar masses mo(l, 2 ) 7 ^ mn(3)fIE\. 


7 Upper bounds from the 19 parameter pMSSM 

The mass bounds from the previous section depend on (reasonable) assumptions about high 
scale physics which are implicit in the NUHM2 model. The previous NUHM2 bounds were 
derived from requiring that no contributions to the renormalization-group-improved one-loop 
scalar potential be far larger than m|/ 2 . The RG-improved scalar potential contains in fact 
leading two loop terms since the parameters entering the scalar potential contain the effects of 
RG running. For instance, the bound on the gluino mass arises mainly from its RG contribution 
to the stop masses: a large value of M 3 pushes m^^ ^ large values leading to large S“(ti 2 ) 
contributions to m^. Likewise, other soft terms contribute to the evolution of to small 
negative values at the weak scale. 

It is popular in recent years to dispense with RG running and examine physics within the 
phenomenological MSSM, a model with 19 free weak scale parameters [36]. Results from the 
pMSSM will look in some cases very different from those obtained with high scale models even 
though in some sense the pMSSM contains models like NUHM 2 , and even though a given 
spectra generated within either of the NUHM 2 or the pMSSM models will yield up exactly 
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the same value of /S.ew- The main difference is that the two-loop contribution to the scalar 
potential arising from RG running of M 3 will no longer occur, thus obviating any bound on 
the gluino mass. In this case, we add in the explicit leading two-loop terms of order atUa as 
computed by Dedes and Slavich[2Z!- These terms contain sensitivity to rrig and to 2 - 

To obtain mass bounds from the pMSSM, we implement a scan over the 19-dimensional 
weak scale parameter space with the following limits: 

• mfR ■ 0.2 - 20 TeV, 

• rui^, rUe^ : 0.1 - 20 TeV, 

• At, Ab, Ar : -40 ^ -h40 TeV, 

• Ml : 0.05 - 10 TeV, M 2 : 0.1 - 10 TeV, M 3 : 0.4 - 10 TeV, 

• p : 100 - 500 GeV, ttia : 0.15 - 20 TeV, 

• tan /3 : 3 — 60. 

We actually input the pMSSM parameters at a scale ApMSSM = 20 TeV- just above the 
maximal soft masses. This ensures the reduced scale dependence of the scalar potential from 
the interplay between the RG running of the soft terms and the log(Q^) dependence of the loop 
contributions S“(z) and The contributions of to Aew are suppressed by tan^/3 

and so typically are of little consequence. 

The pMSSM tree level bounds from Aew < 30 can be read off directly from Eq. [^without 
the need for any scan: 

• /i - 350 GeV, 

• {weak) ~ 350 GeV, 

• ~ 350 GeV tan /5. 

For large then — wia so the latter bound translates to a bound on niA- 

In Fig. we plot the value of Aew vs. the same parameters as in Fig. i a) m-g, b) fi, 
c) Ml and d) M 2 . From frame a) we see a very important result: in the pMSSM, the bound 
on rrig for Aew < 30 has moved up to about 7 TeV, well beyond what was obtained for the 
NUHM2 model. This bound arises after the inclusion of the order atag two-loop contribution 
to the scalar potential which includes sensitivity to rrig |H7|. 

From frame b), we see that, as expected, [i is once again bounded by p ~ 350 GeV for 
Aew < 30, as is required by the tree-level contribution to Eq. In frame c), we plot Aew 
vs. the bino mass Mi. For the pMSSM, mass bounds on Mi arise from the Here, 

we see that Mi can range as high as 9 TeV- far beyond the bounds obtained in the NUHM 2 
model where gaugino mass unification is requred. In that case, the upper limit on rrig ^ 4 TeV 
translates to a bound on Mi ~ 0.8 TeV. Likewise, in frame d) we plot Aew vs. wino mass M 2 . 
In this case, we find M 2 ~ 6 TeV. This bound mainly arises from the contributions S“(ITj). It 
is also far beyond the mass bound from NUHM 2 where M 2 ~ 1.6 TeV was found. 
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Figure 5: Plot of Aew vs. nig, /i, Mi and M 2 from a scan over the 19 weak scale pMSSM 
parameters. Points with Aew < 30 (below dotted line) are considered natural. 
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Figure 6 : Plot of Aew vs. ^62 ^ 19 weak scale pMSSM 

parameters. 


In Fig. [^we display values of vs. a) 171 ^^ h) c) and d) From frame a), 

we find that 7X1^ ^ 3.5 TeV, just slightly larger than the bound arising in the NUHM2 model. 
This bound arises mainly due to the contributions to the scalar potential. In frame 6), 

we find that 10 TeV. This again is slightly beyond the bound arising from the NUHM2 

case with unified matter scalars at ttigut- From frames c) and d), we find that 9 TeV 

and ^10 TeV. These results are in rough accord with values obtained from the NUHM2 
model. 

In Fig. we attempt to extract upper mass bounds on squarks, sleptons, staus and heavy 


Higgs bosons in analogy to Fig. 


From frame a), we find naively that rriu^ ^ 10 TeV. This mass 


bound arises from D-term contributions to first/second generation scalar masses that enter the 
scalar potential via As pointed out in Ref. [71], these D-term contributions all cancel 

amongst themselves provided that one of several mass degeneraey patterns exist 1. separately 
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Figure 7: Plot of ^ew vs. ruf^ and tua from a scan over 19 weak scale pMSSM 

parameters. 

squark and slepton mass degeneracy, 2. separately left- and right- sfermion degeneracy, 3. 
degeneracy within SU (5) multiplets and 4. degeneracy within an entire generation, as expected 
in 5'O(10) GUTs. In these cases, the contributions- which are all proportional to weak isospin 
and hypercharge assignments- necessarily sum to zero for degenerate masses. When these 
contributions thus sum to zero, then there are no bounds on first/second generation squark and 
slepton masses. However, since it is highly improbable to generate these degeneracy patterns 
from a random pMSSM scan, then mass bounds do arise from our random scan. As seen 
from frame a), we expect rriu^ ~ 10 TeV. But this mass bound would disappear if we invoked 
degeneracy conditions amongst the physical masses within any of the patterns listed above. 
Likewise, in frame 6), we see that ~ 10 TeV and from frame c) we expect nif^ ~ 9 TeV. 
The mass bound on tua which arises from the loop contributions E“(h, i7, H^) are shown in 
frame d). These constraints are actually somewhat stronger than the naive tree-level constraint 
of niA ~ 350 GeV tan (3 which is ~ 18 TeV for large tan /5 ~ 50. Instead, we find tua ~ 10 TeV 
from our scan over 19 weak scale pMSSM parameters. 
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Figure 8: Plot of ^ew vs . rrig, rrif^, rrif^ and from a scan over the SUGRA19 model. 

7.1 Comparison between pMSSM and SUGRA19 model 

It may be worthwhile to compare the preceding results from the 19 parameter pMSSM model- 
with inputs defined at the weak scale- with results from Ref. ra where naturalness was 
examined in the context of the SUGRA19 model with 19 input parameters defined at the GUT 
scale. The goal of Ref. ua was to see how how low in ^ew one might go within the context of 
a general SUGRA model, but not to establish upper bounds (which require a thorough rather 
than a focussed parameter space scan). To compare sparticle mass upper bounds between 
pMSSM and SUGRA19, we have re-rnn SUGRA19 inclnding the aforementioned two-loop 
corrections with a thorongh parameter space scan. In both the pMSSM model and the 
SUGRA19 model, the SUSY mu parameter is bounded as p < 350 GeV for ^ew < 30 since 
this quantity enters the scalar potential at tree-level. 

The results from the SUGRA19 model are shown in Fig. [^for /S.ew versus a) nig, h) nif^, 
c) 771^2 and d) . The npper bound on the gluino mass is very similar in the two cases where 
we find rUg ~ 7 TeV. Also, for both the SUGRA19 and pMSSM models, the third generation 
squark masses have similar upper bounds: r\j 3 - 3.5 TeV and - 8 - 10 TeV. 
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From this comparison, and in the spirit of providing a target for experimenters seeking to 
verify or disprove SUSY, we would conclude that a negative search for light higgsinos with mass 
~ 350 GeV would rule out naturalness in the context of the MSSM (more complicated beyond- 
the-MSSM extensions could always be built to circumvent these bounds [23 123 [2S])- Thus, 
non-observation of light higgsinos would rule out both high scale and weak scale renditions of 
SUSY based on the MSSM. Alternatively, failure to find a gluino with mass rrig ~ 4 TeV would 
rule out constrained SUSY models such as NUHM2 with unified soft parameters defined as 
high as Q = mcuT- Colliders with the ability to probe beyond rrig ~ 7 TeV would be needed 
to test/exclude gluino pair production within the pMSSM or SUGRA19 context [HU], 


8 Conclusions: 

In this paper, our goal was to sharpen the upper bounds on sparticle masses arising from 
naturalness in order to provide a target for experimenters seeking to confirm or rule out the 
weak scale supersymmetry hypothesis. 

While most sparticle search results are presented as lower bounds in simplified or complete 
model parameter space, the question arises: how far out in parameter space ought one to go 
before discovering weak scale SUSY or claiming it is dead? While earlier papers by others 
have presented naturalness as a subjective, fuzzy and model-dependent notion, instead here we 
argue that naturalness is 

• objective, 

• model-independent in that different models giving rise to the same spectra have the same 
value of naturalness, and 

• predictive. 

The previous confusion on this subject arose from what constitutes independent model param¬ 
eters. In the case of gravity mediation, for any given hidden sector the soft terms are calculable 
as multiples of the gravitino mass 771 . 3 / 2 , he. they are dependent By appropriately combining 
dependent terms, then the BG measure implies the same general consequences as the model 
independent electroweak measure Aew- We show visually that fine-tuning already arises at 
Aew ~ 20 — 30. To be conservative, we take Aew < 30 to derive upper bounds on parameters 
and sparticle masses. 

In Sec. 1^ we sharpened up previous bounds on sparticle masses by increasing the range 
of parameters enough to ensure that upper bounds arose from Aew < 30 and not from some 
artificial cut-off imposed on our scan limits. For the NUHM2 model which allows SUGRA grand 
unification and values of Aew below 10, we found ~ 4 TeV, well beyond the reach of even 
high luminosity LHC. This bound is much higher than previous estimates which assumed that 
the various soft SUSY breaking terms are independent. Alternatively, we find the superpotential 
mu term to be /i 350 GeV. The range of NUHM2 input parameters which are allowed is shown 
in Fig. 1^ 

The allowed ranges of sparticle masses from NUHM2 are shown in Fig. as the colored 
histograms. As remarked before, the all-important gluino mass can range from nig ~ 1.3 — 4 
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Figure 9: Range of NUHM2 model parameters allowed by naturalness with IS.EW < 30. 

TeV. Third generation squarks, which were argued to exist in the sub-TeV regime in previous 
natural SUSY papers, are found to be consistent with naturalness when mq ~ 3 TeV while 
rrii^ ~ 9 TeV and r\j 10 TeV. A large At parameter here acts to reduce EW fine-tuning 
whilst lifting mu up to ~ 125 GeV. The manifestation of large At is a splitting of the top-squark 
mass eigenstates. First and second generation sfermions can have masses ranging into the 10 
TeV regime so long as they are sufficiently degenerate that naturalness contributions from D- 
terms largely cancel amongst themselves[ZT]. Such heavy matter scalars provide a decoupling 
solution to the SUSY flavor and CP problems and if their mass is comparable to m 3 / 2 , then 
one expects a decoupling solution to the gravitino problem as well. The heavy Higgs bosons 
A, H and can range up to the 5-8 TeV level. While such heavy sparticle and Higgs masses 
are consistent with naturalness, typically we do not expect large deviations from SM rates 
in rare decay branching fraction measurements such as Bg jjA ii~ ov b ^ sj from natural 
SUSY[2n] due to the presence of heavy mediators. 

We have also extracted mass bounds on sparticles by requiring Aew < 30 in a scan over 
the 19 dimensional weak scale parameter set of the pMSSM. These bounds are shown as black 
lines in Fig. [^or as arrows when no bound arises. The pMSSM bounds on scalar masses tend 
to be comparable to those from the NUHM2 model. However, bounds on gaugino masses are 
severely different. By including leading two-loop contributions to the scalar potential, we find 
a bound of rtig < 7 TeV arises in the pMSSM. Also, the bounds on the bino mass Mi ~ 9 TeV 
and wino mass M 2 ~ 6 TeV- which arise from neutralino and chargino loops- are much higher 
than the corresponding upper bounds extracted from the NUHM2 model. A comparison of 
upper bounds extracted from NUHM2 and pMSSM is listed in Table 

Our results have important implications for future particle physics facilities. Even the high 
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Figure 10: Range of sparticle masses allowed by naturalness within NUHM2 model with IS.ew < 
30. The black bars show upper bounds from the pMSSM model with 19 weak scale parameters. 


luminosity LHC can explore only about half of natural SUSY parameter space. However, if we 
require a more stringent naturalness condition of ^ew ~ 10 then rrig ~ 2 TeV and the most 
natural region of parameter space should be accessible to LHC13 (LHC13 has a projected 5cr 
discovery reach to rrig ~ 2 TeV for 300-1000 fb“^ of datajH^). 

The key feature of naturalness- that quasi-degenerate higgsinos lie in the 100-350 GeV mass 
range- highly motivates the construction of an e+e" collider which can operate with ^/s > 2fi. 
Such a machine, constructed initially as a Higgs factory, would turn out to be also a higgsino 
factory which would usher in the era of SUSY discovery while simultaneously elucidating the 
nature of dark matter. In this case, we would expect it to consist of an admixture of higgsino-like 
WIMPs and axionsjnZj. 
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9 Appendix 

In this appendix, we present some details about the radiative contributions E“(z) to the natu¬ 
ralness measure Aew- 
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mass 

NUHM2 

pMSSM 


0.35 

0.35 

(weak) 

0.35 

0.35 

rriA 

5-8 

10 

Ml 

0.9 

9 

M 2 

1.6 

6 

m-g 

4 

7 


3 

3 


9 

9 


9 

9 


10 

10 

nig 

10 (20) 

10/none 

mi 

10 (20) 

10/none 


Table 2: Upper bounds on masses (in TeV) from naturalness with ^ew < 30 from a scan 
over NUHM2 model versus a scan over the 19 weak scale parameter pMSSM. The entries in 
parentheses would result if one allows for non-degenerate generations of soft scalar masses 
mo(l, 2) 7 ^ mo(3)[2S]- The lack of bounds after the slash symbol arise in the case where highly 
degenerate sqnark and slepton masses develop [71] . 


9.1 Limits from E“(ti 2 ) 


Usually, the dominant contributions to come from the top squarks, owing to the large value 
of the top-squark Yukawa coupling ft. The top squark contributions are given by[ 




tl,2 ^ 


ft-glT 


- %l(i - lxw)^t 

m~ — m~ 

t2 


( 10 ) 


where Aj = — m?^)/2 + M| cos 2f3{\ — \xw), xw = sin^ 9w and where 


tR 


m 


F{m ) = m l^log ^ “ ly 

with the optimized scale choice In the denominator of Eq 

expressions for ^ should be used. 


10 


( 11 ) 


the tree level 


In Fig. 


11 


we plot out the S“(ti 2 ) contributions to ^ew in the weak scale vs. At 
plane where we take mi^ = 2.6mf^ which is typical for RNS models. We also adopt tan /3 = 10, 
/i = 150 GeV and ft = 0.8365. In frame a), the blue shaded region bounded by the green 
contour has AEwih) < 30. For At ^ 0, then mi^ ~ 2 TeV although for this value of At it is 
essentially impossible to generate nih as high as 125 GeV[71j. For large stop mixing, i.e. large 
|Af|, then there still exist bands of low AEw{ii) at large \At\ and mi^ ~ 1 TeV which occur 
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AEw{stoPi) AEw(stop2) 



Ai (TeV) At (TeV) 

Figure 11 : Contour plot of a) /S.Ewiii) from and h) AEwih) from ^“(£ 2 ) in the weak- 

scale vs. At plane. The curves correspond to fixed values of AEw(ti, 2 ) as labelled. 


where the mass eigenvalue becomes very small, ~ 100 GeV, so that F(rri~^) is suppressed. 
Alternatively, the low Aevv(^i) bands at large and large At occur due to cancellations in 
the square bracket of Eq. m 

In frame b), where AEwih) is shown, we again see large regions of suppressed fine-tuning 
contributions. Here, mi^ ~ 6 TeV is required for AEw{i 2 ) < 30 although side-bands extend 
out to large \At\ at mi^ ~ 1 and 2.5 TeV. The upper band occurs where mi^/mt^ ~ e so that 
the log term in F(m^^ cancels. The main point is that regions exist with large enough \At\ 
(so that nih —7-125 GeV) and ~ 1 — 3 TeV where both AEwiii, 2 ) become small [2Sl EE]- 
For higher values of one or the other of Aevv(H, 2 ) necessarily becomes > 30, leading to 
fine-tuning. 

Since both S(ti) and ^(^ 2 ) must be small, we show in Fig. 12 the resulting green contour 
from requiring max[AEwiii) i ^Ew{i 2 )] < 30. These combined results show that mi^ is bounded 
from above by about 4.5 TeV. The region below the thick black dashed contour is where 

mi^ < m-t. 


9.2 Limits from 2 ) 


In Fig. we plot the contributions to max[AEwibi), AEwib 2 )] from the 61^2 squarks. We plot 
for tan/5 = 10 and = 0.72m^^ but in the vs. fi plane. For 6 -squark contributions, we 
have 

3 


KM = ,) 


^ 1 , 2 ^ 


2 /bV - - lxw)Ab 

ml — ml 

02 Ol 


( 12 ) 
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Figure 12: Contour plot of max[AEwiii), ^Ewih)] in the weak-scale vs. At plane. The 
region below the thick dashed black contour is where mi^ < nit. 


where A;, = {rn^^ — m?^)/2 — M|cos2/5(| — |xw)- We take fi, — 0.13 and — mi^rrit^ with 


mi^ = 1275 GeV and = 3690 GeV. 


From Fig. 


13 


the contribution A(6i) can become small even for very large ~ 20 TeV, 

but only for large/fine-tuned values of fi. In contrast, AEw(p 2 ) < 30 only for ~ 7 TeV. 
The latter provides a solid upper limit on which is usually ~ 


9.3 Limits from T1^{W\^2) 

The chargino contributions to A^vv are given by 




It 




A 


^ — 2m^ cos 2/3 




W2 


— 


(13) 


Wi 


14 


n 


a), b) in the M 2 vs. fi plane. The 
'act, the AEwiWi) contribution is 


We plot the kFi and IF 2 contributions to Aew in Fig’s 
contribution Aew{^i) is always small for /r ~ 1 TeV. 
small all over the plane except when /i ~ M 2 in which case the denominator in Eq. becomes 
small so that AEwiWi) blows up. In contrast, from AEwiW 2 ) we find that /r, M 2 ~ 5 TeV. 
This illustrates that wino masses all by themselves cannot become too large. 


9.4 Limits from E“(Zi_ 4 ) 

The contributions to Aew from the neutralino mass squared matrix are found to be [20]: 
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Figure 13: Contour plot of AEw{bi,b 2 ) from E“( 6 i^ 2 ) in the weak-scale vs. fi plane. We 
take mi^ = 0.72m^^ and tan/5 = 10. 



Figure 14: Contour plot of AEwiW^i^'i) froni E“(lFi,lF 2 ) in the weak-scale M 2 vs. /r plane 
for tan /5 = 10 . 
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Figure 15: Contour plots of Aew{Zi, Z 2 , Z 3 , Z 4 ) from in the weak-scale M 2 vs. Mi 

plane for fi = 150 GeV and tan (3 = 10. 


where 


K{z^) 


1 

167r2 D{Zi) 


K{Z{) - 2(/ + cos^ (3{m\ 



(14) 


K{Z,) = 




+m; 


-m; 


g\Mt + + g'\Mi + /r^) + {g^ + g'^)M, 


gi\g^Ml + g'^M‘i) + 


/ 2 ^ 


( 15 ) 


D{Zi) — nj^j(m|: ) and m;y — Mi cos^ 6 W+M 2 sin^ 6 w The contribution to max[AEw{Zi)\ 


(z = 1 — 4) are shown in Fig. 15 in the Mi vs. M 2 plane for p = 150 GeV and tan /5 = 10. 
From the hgure we are able to extract that Mi ~ 10 TeV while M 2 ~ 6 TeV. These bounds are 
independent of any assumptions about gaugino mass unification, unlike those of Sec. 


9.5 Limits from 

For Higgs bosons, it is found that [26] 


YZih, H) 


9 % jpt 2 M| + m^(l + 4cos2/5 + 2cos^ 2/5) 


while 
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Figure 16: Contour plot of max[/S.Ew{h^ H, i?^)] from H, H^) in the weak-scale tan P vs. 
niA plane. 


K(H^) = ( 17 ) 

The contributions from each of h, H and are plotted in Fig. [^as max[AEw{h), Aew{H), Aew{H^)] 
in the vs. tan/5 plane. The contribution AEw{h, H, H^) < 30 requires < 8 TeV. The 
largest term usually comes from the charged Higgs contribution. These values are somewhat 
higher than that which comes from the term in Eq. 


References 

[1] H. Baer and X. Tata, “Weak scale supersymmetry: From superfields to scattering events,” Cam¬ 
bridge, UK: Univ. Pr. (2006) 537 p. 

[2] E. Witten, Nucl. Phys. B 188, 513 (1981); R. K. Kaul, Phys. Lett. B 109, 19 (1982). 

[3] G. Aad et al [ATLAS Collaboration], Phys. Lett. B 716 (2012) 1. 

[4] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716 (2012) 30. 

[5] U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B 260, 447 (1991); J. R. Ellis, S. Kelley 
and D. V. Nanopoulos, Phys. Lett. B 260 (1991) 131; P. Langacker and M. x. Luo, Phys. Rev. 
D 44 (1991) 817. 

[6] L. E. Ibanez and G. G. Ross, Phys. Lett. BllO, 215 (1982); K. Inoue et al. Prog. Theor. Phys. 68, 
927 (1982) and 71, 413 (1984); L. Ibanez, Phys. Lett. B118, 73 (1982); H. P. Nilles, M. Srednicki 


27 










and D. Wyler, Phys. Lett. B 120 (1983) 346; J. Ellis, J. Hagelin, D. Nanopoulos and M. Tamvakis, 
Phys. Lett. B125, 275 (1983); L. Alvarez-Gaume. J. Polchinski and M. Wise, Nucl. Phys. B221, 
495 (1983); B. A. Ovrut and S. Raby, Phys. Lett. B 130 (1983) 277; for a review, see L. E. Ibanez 
and G. G. Ross, Comptes Rendus Physique 8 (2007) 1013. 

[7] M. S. Carena and H. E. Haber, Prog. Part. Nucl. Phys.bO (2003) 63 [hep-ph/0208209]. 

[8] G. Aad et al. (ATLAS collaboration), Phys. Lett. B 710 (2012) 67. 

[9] S. Chatrchyan et al. (CMS collaboration), Phys. Rev. Lett. 107 (2011) 221804. 

[10] R. Barbieri and G. Giudice, Nucl. Phys. B 306 (1988) 63. 

[11] G. W. Anderson and D. J. Castano, Phys. Rev. D 52 (1995) 1693. 

[12] J. R. Ellis, K. Enqvist, D. V. Nanopoulos and F. Zwirner, Mod. Phys. Lett. A 1 (1986) 57. 

[13] B. C. Allanach, K. Cranmer, C. G. Lester and A. M. Weber, JHEP 0708 (2007) 023; M. E. Cabr¬ 
era, J. A. Casas and R. Ruiz de Austri, JHEP 0903 (2009) 075; M. E. Cabrera, J. A. Casas and 
R. Ruiz d Austri, JHEP 1005 (2010) 043; S. Fichet, Phys. Rev. D 86 (2012) 125029; D. Kim, 
P. Athron, C. Balzs, B. Farmer and E. Hutchison, Phys. Rev. D 90 (2014) 5, 055008; A. Fowlie, 
Phys. Rev. D 90 (2014) 1, 015010; A. Fowlie, arXiv: 1506.03786 [hep-ph]. 

[14] S. S. AbdusSalam, Int. J. Mod. Phys. A 29 (2014) 27, 1450160. 

[15] S. S. AbdusSalam and L. Velasco-Sevilla, arXiv: 1506.02499 [hep-ph]. 

[16] Joint LEP 2 Supersymmetry Working Group, Combined LEP Chargino Results up to 208 GeV, 
http: //lepsusy. web. cern. ch/lepsusy/www/inosnnoriond01/charginos_pub. html. 

[17] R. Barbieri and A. Strumia, Phys. Lett. B 433 (1998) 63. 

[18] J. Lykken and M. Spiropulu, Sci. Am. 310N5 (2014) 5, 36. 

[19] S. Dimopoulos and L. Susskind, Nucl. Phys. B 155 (1979) 237. 

[20] S. Weinberg in To Explain the World: The Discovery of Modern Science (Harper, 2015, 416p). 

[21] H. Baer, V. Barger, D. Mickelson and M. Padeffke-Kirkland, Phys. Rev. D 89 (2014) 115019. 

[22] K. L. Chan, U. Chattopadhyay and P. Nath, Phys. Rev. D 58 (1998) 096004; S. Akula, M. Liu, 
P. Nath and G. Peim, Phys. Lett. B 709 (2012) 192; M. Liu and P. Nath, Phys. Rev. D 87 (2013) 
9, 095012. 

[23] M. Dine, A. Kagan and S. Samuel, Phys. Lett. B 243 (1990) 250; A. Cohen, D. B. Kaplan and 
A. Nelson, Phys. Lett. B 388 (1996) 588; N. Arkani-Hamed and H. Murayama, Phys. Rev. D 56 
(1997) R6733; T. Moroi and M. Nagai, Phys. Lett. B 723 (2013) 107. 

[24] S. Weinberg, Phys. Rev. Lett. 48 (1982) 1303; M. Y. Khlopov and A. D. Linde, Phys. Lett. B 
138 (1984) 265; M. Kawasaki, K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 78 (2008) 
065011. 


28 


[25] H. Baer, V. Barger, P. Huang, A. Mustafayev and X. Tata, Phys. Rev. Lett. 109 (2012) 161802. 

[26] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev and X. Tata, Phys. Rev. D 87 (2013) 
11, 115028. 

[27] P. J. Fox, A. E. Nelson and N. Weiner, JHEP 0208 (2002) 035; A. E. Nelson and T. S. Roy, 
Phys. Rev. Lett. 114 (2015) 201802. 

[28] T. Cohen, J. Kearney and M. Luty, Phys. Rev. D 91 (2015) 075004. 

[29] S. P. Martin, Phys. Rev. D 92 (2015) 3, 035004. 

[30] N. Arkani-Hamed, T. Han, M. Mangano and L. T. Wang, arXiv:1511.06495 [hep-ph]. 

[31] J. E. Kim and H. P. Nilles, Phys. Lett. B 138 (1984) 150. 

[32] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B 104, 199 (1981); A. P. Zhitnitskii, Sov. J. 
Phys. 31, 260 (1980). 

[33] E. J. Chun, Phys. Rev. D 84, 043509 (2011); K. J. Bae, E. J. Chun and S. H. Im, JCAP 1203, 
013 (2012). 

[34] H. Murayama, H. Suzuki and T. Yanagida, Phys. Lett. B 291, 418 (1992); 

[35] K. J. Bae, H. Baer and H. Serce, Phys. Rev. D 91, 015003 (2015). 

[36] A. Djouadi, J. -L. Kneur and G. Moultaka, Comput. Phys. Commun. 176 (2007) 426; J. Conley, 
S. Gainer, J. Hewett, M. Le and T. Rizzo, Eur. Phys. J. C 71 (2011) 1697; S. S. AbdusSalam, 
AIP Conf. Proc. 1078 (2009) 297; S. S. AbdusSalam, B. C. Allanach, F. Quevedo, F. Feroz 
and M. Hobson, Phys. Rev. D 81 (2010) 095012; S. Sekmen, S. Kraml, J. Lykken, F. Moortgat, 
S. Padhi, L. Pape, M. Pierini and H. B. Prosper et al.^ J. High Energy Phys. 1202 (2012) 075; 
M. W. Cahill-Rowley, J. L. Hewett, A. Ismail and T. G. Rizzo, Phys. Rev. D 86 (2012) 075015 
and Phys. Rev. D 88 (2013) 3, 035002; C. Boehm, P. S. B. Dev, A. Mazumdar and E. Pukartas, 
JHEP 1306 (2013) 113; A. Arbey, M. Battaglia, A. Djouadi and F. Mahmoudi, Phys. Lett. B 
720 (2013) 153; M. Chakraborti, U. Chattopadhyay, A. Choudhury, A. Datta and S. Poddar, 
JHEP 1407 (2014) 019. 

[37] A. Dedes and P. Slavich, Nucl. Phys. B 657 (2003) 333. 

[38] A. Mustafayev and X. Tata, Indian J. Phys. 88 (2014) 991 jarXiv: 1404.1386 [hep-ph]]. 

[39] L. E. Ibanez, C. Lopez and C. Munoz, Nucl. Phys. B 256 (1985) 218; A. Lleyda and C. Munoz, 
Phys. Lett. B 317 (1993) 82. 

[40] H. Abe, T. Kobayashi and Y. Omura, Phys. Rev. D 76 (2007) 015002. 

[41] S. P. Martin, Phys. Rev. D 75 (2007) 115005. 

[42] J. L. Feng, Ann. Rev. Nucl. Part. Sci. 63 (2013) 351. 

[43] S. K. Soni and H. A. Weldon, Phys. Lett. B 126 (1983) 215. 


29 


[44] V. S. Kaplunovsky and J. Louis, Phys. Lett. B 306 (1993) 269. 

[45] A. Brignole, L. E. Ibanez and C. Munoz, Nucl. Phys. B 422 (1994) 125 [Erratum-ibid. B 436 
(1995) 747]. 

[46] H. Baer, V. Barger and M. Savoy, Phys. Scripta 90 (2015) 6, 068003. 

[47] R. Harnik, G. D. Kribs, D. T. Larson and H. Murayama, Phys. Rev. D 70 (2004) 015002. 

[48] R. Kitano and Y. Nomura, Phys. Lett. B 631 (2005) 58 and Phys. Rev. D 73 (2006) 095004. 

[49] M. Papucci, J. T. Ruderman and A. Weiler, JHEP 1209 (2012) 035; C. Brust, A. Katz, 
S. Lawrence and R. Sundrum, JHEP 1203 (2012) 103; J. A. Evans, Y. Kats, D. Shih and 

M. J. Strassler, JHEP 1407 (2014) 101. 

[50] S. Cassel, D. M. Ghilencea, S. Kraml, A. Lessa and G. G. Ross, JHEP 1105 (2011) 120; 
G. G. Ross, Eur. Phys. J. C 74 (2014) 2699. 

[51] H. P. Nilles, Phys. Lett. B 115 (1982) 193; A. Chamseddine, R. Arnowitt and P. Nath, Phys. 
Rev. Lett. 49 (1982) 970; R. Barbieri, S. Ferrara and C. Savoy, Phys. Lett. B119 (1982) 343; 

N. Ohta, Prog. Theor. Phys. 70, 542 (1983); L. Hall, J. Lykken and S. Weinberg, Phys. Rev. 
D27 (1983) 2359; for an early review, see e.g. H. P. Nilles, Phys. Rept. 110 (1984) 1; for a recent 
review, see D. J. H. Chung, L. L. Everett, G. L. Kane, S. F. King, J. D. Lykken and L. T. Wang, 
Phys. Rept. 407 (2005) 1. 

[52] G. Kane, C. Kolda, L. Roszkowski and J. Wells, Phys. Rev. D 49 (1994) 6173. 

[53] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev and X. Tata, Phys. Rev. D 87 (2013) 
3, 035017. 

[54] H. Baer, V. Barger and D. Mickelson, Phys. Rev. D 88 (2013) 9, 095013. 

[55] R. D. Peccei, Lect. Notes Phys. 741 (2008) 3. 

[56] L. F. Abbott and P. Sikivie, Phys. Lett. B120 (1983) 133; J. Preskill, M. Wise and F. Wilczek, 
Phys. Lett. B120 (1983) 127; M. Dine and W. Fischler, Phys. Lett. B120 (1983) 137; M. Turner, 
Phys. Rev. D33 (1986) 889. 

[57] K. J. Bae, H. Baer and E. J. Chun, Phys. Rev. D 89 (2014) 3, 031701; K. J. Bae, H. Baer and 
E. J. Chun, JCAP 1312 (2013) 028; K. J. Bae, H. Baer, V. Barger, M. R. Savoy and H. Serce, 
Symmetry 7 (2015) 2, 788 jarXiv: 1503.04137 [hep-ph]]. 

[58] H. Baer, V. Barger, P. Huang, D. Mickelson, M. Padeffke-Kirkland and X. Tata, Phys. Rev. D 
91 (2015) 7, 075005. 

[59] ISAJET 7.84, by H. Baer, F. Paige, S. Protopopescu and X. Tata, hep-ph/0312045. 

[60] H. Baer, C. H. Chen, R. Munroe, F. Paige and X. Tata, Phys. Rev. D 51 (1995) 1046; H. Baer, 
J. Ferrandis, S. Kraml and W. Porod, Phys. Rev. D 73 (2006) 015010. 


30 


[61] D. Matalliotakis and H. P. Nilles, Nucl. Phys. B 435 (1995) 115; P. Nath and R. L. Arnowitt, 
Phys. Rev. D 56 (1997) 2820; J. Ellis, K. Olive and Y. Santoso, Phys. Lett. B 539 (2002) 107; J. 
Ellis, T. Falk, K. Olive and Y. Santoso, Nucl. Phys. B 652 (2003) 259; H. Baer, A. Mustafayev, 
S. Profumo, A. Belyaev and X. Tata, J. High Energy Phys. 0507 (2005) 065. 

[62] H. Baer, V. Barger, A. Lessa and X. Tata, Phys. Rev. D 86 (2012) 117701 jarXiv: 1207.4846 
[hep-ph]]. 

[63] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev, W. Sreethawong and X. Tata, JHEP 
1312 (2013) 013. 

[64] B. Altunkaynak, H. Baer, V. Barger and P. Huang, arXiv: 1507.00062 [hep-ph]. 

[65] H. Baer, V. Barger and P. Huang, JHEP 1111 (2011) 031. 

[66] Z. Han, G. D. Kribs, A. Martin and A. Menon, Phys. Rev. D 89 (2014) 7, 075007. 

[67] H. Baer, A. Mustafayev and X. Tata, Phys. Rev. D 90 (2014) 11, 115007. 

[68] C. Han, D. Kim, S. Munir and M. Park, JHEP 1504 (2015) 132. 

[69] H. Baer, V. Barger, D. Mickelson, A. Mustafayev and X. Tata, JHEP 1406 (2014) 172. 

[70] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev, W. Sreethawong and X. Tata, Phys. 
Rev. Lett. 110 (2013) 15, 151801. 

[71] H. Baer, V. Barger, M. Padeffke-Kirkland and X. Tata, Phys. Rev. D 89 (2014) 3, 037701. 

[72] H. Baer, V. Barger and M. Padeffke-Kirkland, Phys. Rev. D 88 (2013) 055026. 

[73] H. Baer, C. Balazs, P. Mercadante, X. Tata and Y. Wang, Phys. Rev. D 63 (2001) 015011. 

[74] H. Baer, V. Barger and A. Mustafayev, Phys. Rev. D 85 (2012) 075010. 


31 


